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We report the measurement of the Cabibbo-Kobayashi-Maskawa CP-violating angle 7 through a 
Dahtz plot analysis of neutral D meson decays to K^-k'^-k^ and K'^K^ K~ produced in the processes 
B^ -^ DK^, B^ -^ D*K^ with D* -^ Dtt", D-/, and B"^ -^ DK'*^ with K*^ -^ KI-k^ , using 468 
million BB pairs collected by the BABMl detector at the PEP-II asymmetric-energy e^e^ collider 
at SLAC. We measure 7 = (68 ± 14 ± 4 ± 3)° (modulo 180°), where the first error is statistical, 
the second is the experimental systematic uncertainty and the third reflects the uncertainty in the 
description of the neutral D decay amplitudes. This result is inconsistent with 7 = (no direct CP 
violation) with a significance of 3.5 standard deviations. 

PACS numbers; 13.25.Hw, ll.30.Er, 12.15.Hh, 13. 25. Ft 

The breaking of the CP symmetry in the quark sec- only tree-level or both tree- and loop-level amplitudes, 

tor of the electroweak interactions arises in the stan- The comparison of these two classes of measurements 

dard model (SM) from a single irreducible phase in the tests the CKM mechanism, thus offering a strategy to 

Cabibbo-Kobayashi-Maskawa (CKM) quark- mixing ma- search for new physics [2]. The angle 7 of the unitar- 

trix [1]. This phase can be measured using a variety of ity triangle, defined as arg [— ^d^*h/V'gdV'ch]' where Vij 

methods involving B- meson decays mediated by either are elements of the CKM matrix, is particularly relevant 



since it is the only CP-violating parameter that can be 
cleanly determined using solely tree-level i3-meson de- 
cays. Its precise determination constitutes an important 
goal of present and future experiments in flavor physics. 

In B^ -^ DK^ decays [3, 4] the color- favored B~ ^i- 
D°K~ (6 -^ cus) and the color-suppressed B~ -^ D^K~ 
{b — >■ ucs) transitions [5] interfere when the D'^ and 
D° decay to a common final state [6]. The two inter- 
fering amplitudes differ by a factor rBe^^^'^^''\ where 
rB is the magnitude of the ratio of the amplitudes 
AiB- -^ WK-) and A{B- -^ D^R-), and 5b is 
their relative strong phase. An amplitude analysis of 
the Dalitz plot (DP) of D" and D^ mesons decaying into 
the iCgTr+TT^ and K^K'^K~ self-conjugate final states 
from B^ -^ DK^ decays offers a unique way to ac- 
cess the complex amplitude ratios and thus the weak and 
strong phases, and tb- The experimental sensitivity to 
7 arises mostly from regions in the DP where Cabibbo- 
favored (CF) and doubly-Cabibbo-suppressed (DCS) am- 
plitudes interfere, and from regions populated by CP 
eigenstates, thus the uncertainty on 7 depends on I/tb 
(rs- 0.1 -0.2). 

In this Letter we study the interference between color- 
favored and color-suppressed transitions as a function 
of the position in the DP of squared invariant masses 



m^K'ih- 



s+ 



'{K^h'^), where h repre- 



sents TT or K, for three related B decays, B^ — > DK^, 
B^ -^ D*K^, and 5=F -^ DK*^ [4, 7], and report the 
most precise single measurement of the complex ampli- 
tude ratios and evidence for direct CP violation. We use 
the complete data sample of 425 fb^^ of integrated lu- 
minosity at the T{4:S), corresponding to 468 x 10^ BB 
pairs, and 45 fb^ at a center-of-mass (cm.) energy 
40 MeV below the T{4S), recorded by the BaBAR exper- 
iment [8] at the PEP-II asymmetric-energy e^e~ collider 
at SLAC from 1999 to 2008. This measurement updates 
our previous results based on a partial sample of 383 x 10® 
BB pairs, from which we reported a significance of direct 
CP violation (7 7^ 0) of 3.0 standard deviations, while 
most of the analysis details remain unchanged [9] . The 
Belle Collaboration using B=F -^ D^*^ K^ , D -)■ K^tt+tt- 
alone [10] has also reported 7 7^ with a significance of 
3.5 standard deviations. 

We reconstruct a total of eight signal samples, B^ -^ 
D^*^K^ and B^ -^ DK*^, with D* -^ DTr°,Dj, 
K*^ — > K^Itt^ , with selection criteria nearly identical 
to our previous analysis. The DK*^ final state, for 
D — > K'^K^K^ , has been considered for the first time. 
For Kg — > TT+TT^ candidates, we further require the de- 
cay length (defined by the K'^ production and decay ver- 
tices) projected along the K'^ momentum to be greater 
than 10 times its error. This additional requirement 
helps to reduce to a negligible level background events 
from D -^ 7r+7r^/i+/i^ decays, and from ai(1260)^ 
misreconstructed as K*^. After all the selection crite- 
ria the background is completely dominated by random 



combinations of tracks arising from continuum events, 
e^e~ -^ qq {q — u,d,s, or c). Background contributions 
from D — > K'^Kg decays are found to be negligible. The 
B^ candidates are characterized using the beam-energy 
substituted B mass rriES, the difference between the re- 
constructed energy of the B^ candidate and the beam 
energy in the e~^e~ cm. frame AE, and a Fisher dis- 
criminant T that combines four topological variables op- 
timized to separate continuum events [9] . We retain can- 
didates with the loose requirements toes > 5-2 GeV/c^, 
-80 < AE < 120 MeV, and \T\ < 1.4, which provide sig- 
nal and sideband regions while removing poorly recon- 
structed candidates [11]. The reconstruction efficien- 
cies in a signal region with toes > 5.272 GeV/c^ and 
|A£;| < 30 MeV are 26%, 12%, 15%, and 14%, for the 
DK'f', D*[DT:^]K'f', D*[D-f]K'f', and DK*'^' final states, 
respectively, for D — > K'^tt^tt^ (and slightly lower for 
D -^ KlK+K-). These values are about 30%, 40%, 
30%, and 20% larger than in our previous analysis, with 
similar background levels, reflecting improvements in 
tracking and particle identification. The to-es, AE, J-, 
and {s-,Sj[.) distributions for events in the signal region 
can be found in [11]. 

The D^ -^ K'^h^h^ decay amplitudes A{s^,s+) are 
determined using the same data sample through DP anal- 
yses of I?'^ mesons from _D*+ — > Z)°7r+ decays produced in 
e+e^ — > cc events [9, 12]. The charge of the low momen- 
tum 7r+ from the £)*+ decay identifies the flavor of the 
D meson. The signal purities of the samples are 98.5% 
and 99.2%, with about 541000 and 80 000 candidates, 
for K'^tt^tt" and K'^K^K~ , respectively. The dynami- 
cal properties of the P- and D-wave amplitudes are pa- 
rameterized through intermediate resonances with mass- 
dependent relativistic Breit-Wigner (BW) or Gounaris- 
Sakurai propagators, Blatt-Weisskopf centrifugal barrier 
factors, and Zemach tensors for the angular distribu- 
tions [13]. The TTTT S-wave dynamics is described through 
a K-matrix formalism with the P-vector approximation 
and 5 poles [9, 14]. For the Ktt S-wave we include a BW 
for the ii'Q(1430)^ state with a coherent non-resonant 
contribution parameterized by a scattering length and 
effective range similar to those used to describe Ktt scat- 
tering data [15]. For the KK S-wave, a coupled-channel 
BW is used for the ao(980) with single BWs for /o(1370) 
and ao(1450) states. Overall, the amplitude models re- 
produce well the DP distributions [12]. MC studies show 
that a significant contribution to the discrepancies arise 
from imperfections modeling the efficiency variations at 
the boundaries of the DP and the invariant mass resolu- 
tion. We account for these and other imperfections in the 
modeling of the D'^ decay amplitudes through our model 
systematic uncertainties. 

We perform a simultaneous, unbinned, and extended 
maximum-likelihood fit (referred to as CP fit) to the 



.(*) 



and 



BT ^ d(*^K^ and B=F ^ DK*^ decay rates F^ 

T s^ as a function of toes, ^E, J-, and (5^,5+) [9, 11] 



We extract the signal and background yields, along with consistent results for all the CP- violating parameters [11] 



the CP-violating parameters z 
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defined as the B^ complex ampli- 



tude ratios z-^' = 
respectively. Here, 
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' BT 



and 



ST — ST 

s^ are the corresponding 
u and h ^f c am- 
the relative strong 



magnitude ratios between the h —, 
plitudes for B^ decays, 5g and (5, 
phases, and k an effective hadronic parameter that ac- 
counts for the interference between B^ -^ DK*^ and 
other B^ — > DK'^tt^ decays, as a consequence of the 
K*^ natural width [9, 16, 17]. Assuming no CP viola- 
tion and neglecting Z)° — 1)° mixing in I?° — > K^h'^h'' 
decays [12, 18, 19], the relation A{s-,s+) = A{s+,s-) 
holds, where A is the D° decay amplitude. The B^ -^ 
D^*^K^ (and similarly for B^ -^ DK*^ replacing z^' 



and r^^ by z 
are then 



sip and rgzp , respectively) signal decay rates 



rW(._ 



,s+) oc 1^^ 



(*)■ 

'b^ 



\A, 



2Az^*U^^i, 



with Azfz = A{szfz,s±), and A = -|-1 except for B^ — > 
D*[Dj]K^ where A = — 1 [20]. We apply corrections for 
efficiency variations and neglect the invariant mass res- 
olution across the DP [9]. For each signal sample, the 
following background components are considered: con- 
tinuum events, B^ — > D^*'tt^ decays where the pion is 
misidentified as a kaon (only for B^ -^ D'^*'^ K^ sam- 
ples), and T{AS) -)■ BB (other than B^ -^ D'^*^tt^) 
decays. The reference CP fit requires events to satisfy 
|Ai?| < 30 MeV, but alternative fits are performed vary- 
ing the requirements on the ttieSj ^E, and T variables 
(e.g. mEs > 5.272 GeV/c^ or J" > -0.1) to study the 
stability of the results. The probability density func- 
tions (PDFs) introduced to describe the signal, contin- 
uum, and K/tt misidcntification components, along with 
the K/iT misidcntification yields, are determined using 
events from signal and B^ — > I?(*V^,I?ai(1260)^ con- 
trol samples. The PDFs for BB background events are 
obtained from large Monte Carlo (MC) samples with full 
detector simulations [9]. 

The CP fit yields 896 ± 35 (154 ± 14), 255 ± 21 
(56 ± 11), 193 ± 19 (30 ± 7), and 163 ± 18 (28 ± 6) signal 
DK^, D*[Dtt^]K^, D*[D"f]K^, and DK*^ events, re- 
spectively, for the X^Tr+TT" [KlK+K^) final state. The 
results for the CP-violating parameters z^!'^ and Zs± are 
summarized in Table I. Figure 1 shows the 39.3% and 
86.5% 2-dimensional confidence-level (CL) contours in 
the z^, z* , and Zgzp planes, corresponding to one- and 
two-standard deviation regions, including statistical er- 
rors only. The distance between the z_ and z+ cen- 
tral values (and similarly for z* and z,,^) is equal to 
2rgT\ sin7|, and the angle defined by the lines connect- 
ing the central values with the origin is 27, and thus is a 
measurement of direct CP violation. Fitting separately 
the data for XStt+tt" and K^K^K^ final states we find 



The first 



TABLE I: CP-violating complex parameters z^ = a;!^ 
and ZsT — Xs^i + iys^ as obtained from the CP fit. 
error is statistical, the second is the experimental systematic 
uncertainty and the third is the systematic uncertainty asso- 
ciated with the D" decay ampHtude models. 





Real part (%) 


Imaginary part (%) 


z_ 


6.0 ±3.9 ±0.7 ±0.6 


6.2 ± 4.5 ±0.4 ±0.6 


z+ 


-10.3 ±3.7 ±0.6 ±0.7 


-2.1± 4.8 ±0.4 ±0.9 


zl 


-10.4 ±5.1 ±1.9 ±0.2 


-5.2± 6.3 ±0.9 ±0.7 


A- 


14.7 ±5.3 ±1.7 ±0.3 


-3.2± 7.7 ±0.8 ±0.6 


Zs- 


7.5 ±9.6 ±2.9 ±0.7 


12.7 ± 9.5 ±2.7 ±0.6 


Zs+ 


-15.1 ±8.3 ±2.9 ±0.6 


4.5 ±10.6 ±3.6 ±0.8 



Experimental systematic errors [9, 11] originate from 
uncertainties in the description of the efficiency varia- 
tions across the DP, the modeling of the DP distribu- 
tions for background events containing misreconstructed 
D mesons, the fractions of continuum and BB back- 
ground events containing a real D meson with either a 
negatively- or positively-charged kaon (or K*), and from 
residual direct CP violation in the B^ — ;> D'^*V± and 
BB background components. We also account for statis- 
tical and systematic uncertainties in the toeSj Ai?, and 
T PDF shapes for signal and background components, 
and the K/n misidcntification yields. These uncertain- 
ties account for effects that arise from the dependence 
of the Toes and F PDF shapes on the chosen AE signal 
region, the differences in BB background for real and 
misreconstructed D mesons, and our limited knowledge 
of the Toes endpoint, the peaking contributions to the 
small BB background, and the e^e~ cm. frame. Smaller 
systematic uncertainties originate from the DP resolu- 
tion, wrongly reconstructed signal events with a real D 
and a kaon (or K*) from the other B meson decay, the 
selection of B candidates sharing tracks with other can- 
didates, and numerical precision in the evaluation of the 
PDF integrals. We also account for residual cross-feed 
of PT ^ D*[D'K°]K'f events into the P± -^ D*[D-f]K'^ 
sample (about 5%), and the estimated uncertainty on the 
hadronic parameter k = 0.9 ±0.1 in the P± — s- DK*^ 
sample [9, 21]. 

Assumptions in the £>" decay amplitude models are 
also a source of systematic uncertainty [9, 11, 12]. We 
use alternative A{s^, s+) models where the BW parame- 
ters are varied according to their uncertainties or within 
the ranges allowed by measurements from other experi- 
ments, the reference K-matrix solution [9] is replaced by 
other solutions [14], and the standard parameterizations 
are substituted by other related choices. These include 
replacing the Gounaris-Sakurai and Kit S-wave parame- 
terizations by BW lineshapes, removing the mass depen- 
dence in the P-vector [22] , changes in form factors such as 
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FIG. 1: (color online). Contours at 39.3% (dark) and 86.5% (light) 2-dimensional CL in the (a) z^, (b) z^, and (c) Zsip planes, 
corresponding to one- and two-standard deviation regions (statistical only), for B~ (solid lines) and B^ (dotted lines) decays. 



changes in the Blatt-Weisskopf radius, and adopting a he- 
licity formalisni [13] to describe the angular dependence. 
Other models are built by removing or adding resonances 
with small or negligible fractions. We find that the over- 
all amplitude model uncertainty on the CP parameters 
are dominated by alternative models built to account for 
experimental systematic effects in the determination of 
A{s-,s+) using tagged D mesons [12]. The statistical 
errors and variations in the A{s^, 5+) model parameters 
with and without D'^-D^ mixing are also propagated to 
z*^ and Zszp. 

Experimental and amplitude model systematic uncer- 
tainties [11] have been reduced with respect to our previ- 
ous measurement [9] as consequence of the use of larger 
data and Monte Carlo samples, and the smaller experi- 
mental systematic contributions to the model uncertainty 
resulting from the improvements in the analysis of tagged 
D mesons [12]. 

A frequentist construction of 1-dimensional confidence 
intervals of the physically relevant parameters p = 
{'^,rB,r*g,Krs,5B,5*B,5s) based on the vector of mea- 
surements z = (z_, z+, zl, z^, Zs_, Zs4.) and their corre- 
lations [11] has been adopted [9]. The procedure takes 
into account unphysical regions which may arise since we 
allow B~ and B^ events to have different r*!^ , Tszp in the 
z measurements. Figure 2 shows 1 — CL as a function of 
7 for each of the three B decay channels separately and 



their combination. Similar scans for r^ , Kr^, Sg , 



^*) 



r(*) 



and 



Ss can be found in [11]. The method has a single ambi- 
guity in the weak and strong phases. The results for all 
the p parameters are listed in Table II. The significances 
of direct CP violation (7 7^ 0) are 1-CL=6.8 x IQ-^, 
5.4 X 10-3, 6.3 X 10"^ and 4.6 x 10"'', which corre- 
spond to 2.7, 2.8, 1.9, and 3.5 standard deviations, for 
B^ -^ DK^, B^ -^ D*K^, B^ -> DK*^, and their 
combination, respectively. 

We have presented a measurement of the b ^ u 
to 6 — > c complex amplitude ratios in the processes 
ST ^ d(*'^K^ and B=F -^ DK*^, using a combined DP 




-150 -100 -50 50 100 150 

Y (deg) 

FIG. 2: (color online). 1 — CL as a function of 7 for 
B^ -> DK"^, B^ ^ D*K'^, and B^ -^ DK*"^ decays sep- 
arately, and their combination, including statistical and sys- 
tematic uncertainties. The dashed (upper) and dotted (lower) 
horizontal lines correspond to the one- and two-standard de- 
viation intervals, respectively. 



TABLE II: The 68.3% and 95.4% 1-dimensional CL regions, 
equivalent to one- and two-standard deviation intervals, for 
7, 5g , 5s, r^ , and K,rs, including all sources of uncertainty. 
The values inside {} brackets indicate the symmetric error 
contributions to the total error coming from experimental and 
amplitude model systematic uncertainties. 



Parameter 


68.3% CL 


95.4% CL 


in 


68^1^4,3} 


[39, 98] 


TB (%) 


9.6 ±2.9 {0.5,0.4} 


[3.7, 15.5] 


r*B (%) 


13.3l3g {1-3,0.3} 


[4.9,21.5] 


Krs (%) 


14.9+^-^ {2.6,0.6} 


<28.0 


5b n 


119t^;^ {3,3} 


[75, 157] 


5b n 


-82 ±21 {5,3} 


[-124,-38] 


5s n 


111 ±32 {11,3} 


[42, 178] 



analysis of D ^' K^tt^tt~ and D — >■ K'^K~^K~ decays. 
The results have improved precision and are consistent 
with our previous measured values [9] and with those 
reported by the Belle Collaboration with D — >■ i^^Tr+Tr" 
alone [10], and with determinations based on other D me- 
son final states [21, 23, 24]. From our measurement we 
determine 7 = (68± 14±4±3)° (modulo 180°), exclude 
the no direct CP- violation hypothesis (i.e., 7 = 0) with 
a CL equivalent to 3.5 standard deviations, and derive 
the most precise single determinations of the magnitude 
ratios r^ and nvg. 
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FIG. 1: (color online). The tties (first column), AE (second column), and T (third column) distributions for (a)-(c) B^ — )■ 
DK^, (d)-(f) B^ -^ D*[L>7r°]7^=F, (g)-(i) B^ ^ D*[D'y]K'^, and (j)-(l) B^ -^ DK'^ decays, with D -f KItv+tv- . The 
distributions are for events in the signal region defined through the requirements jties > 5.272 GeV/c^, |AJ5| < 30 MeV, and 
T > —0.1, except the one on the plotted variable, after all the selection criteria are applied. The curves superimposed represent 
the projections of the CP fit: signal plus background (solid black lines), the continuum plus BB background contributions 
(dotted red lines), and the sum of the continuum, BB, and K/n misidentification background components (dashed blue lines). 
The reconstruction efliciencies (purities) in the signal region, based on simulation studies, are 22% (68%), 10% (81%), 12% 
(55%), and 12% (58%), respectively. 




d) 



5.25 

m^ (GeV/c^) 



5.2 



j,irtBi,ul:irtKu a i.u.i.ij: 




5.25 

mEs(GeV/c') 




0.05 0.1 
AE (GeV) 



1 


g) 


s 




~ 10 


, 


^ 


■'ii 






Ui e 


, A 






1 


1 


w 


„ 


iin 


w 


h llr 


irrr ,tH 



5.2 



5.25 

m,3(GeV/c') 



> lu 




i 


1) 




K? 






a 






c 
























m 










5 - 


■■ 


■ 


■■ 






oO 


in 


, , 1 


ffff 


ll 


•[rr 



-0.05 



0.05 0.1 
AE (GeV) 



§ 

I 10 
UJ 


i) 




5 


- 




:■•' 


J, ... I 


1 






u 


Ift 


w 


.^4 


III II 



1 

Fisher 




5.25 

itIes (GeV/c^) 




-l> 




" n/1 


-■ 






,f 11 


\m 


" 


mil 


-1 




D 


1 
Fi 


sher 



FIG. 2: (color online). Same as in Fig. 1 but for (a)-(c) B^ -^ DK^, (d)-(f) B=F ^ D*[D7V°]K'^, (g)-(i) B^ -> O'lD-ylK"^, 
and (j)-(l) -B^ — >■ DK*^ decays, with D — >■ KgK'^ K~ . The reconstruction efficiencies (purities) in the signal region, based on 
simulation studies, are in this case 20% (82%), 9% (87%), f2% (78%), and If% (81%), respectively. 
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FIG. 3: (color online). The DP distributions for (a)(b) B^ ^ DA'^, (c)(d) B^ ^ D*[D7r°]A:T, (e)(f) B^ -> D*[D7]A'T, 
and (g)(h) B"^ ->■ DK*^ decays, with D -)> KsTt+tv' (left panel) and D ->• AgA'+A'" (right panel). The distributions are 
for events in the signal region defined through the requirements tties > 5.272 GeV/c^, |A_B| < 30 MeV, and T > —0.1, after 
all the selection criteria are applied, and are shown separately for B~ (first and third columns) and B"*" (second and last 
column) decays. For B~ and B^ decays the variables S- and s+ are interchanged. The contours (solid red lines) represent the 
kinematical limits of the D decay. 



TABLE I: Summary of the main contributions to the experimental systematic error on the CP parameters. All contributions 
have been evaluated using the same procedure as in our previous analysis [9]. The statistical contribution to the total error 
has been decreased, as consequence of the use of larger data and Monte Carlo (with full detector simulation) samples. For 
example, larger simulated continuum samples help to significantly reduce the uncertainty arising from the modeling of the DP 
distributions for background events containing misreconstructed D mesons. 

Source x- y- x+ y+ x*_ y*_ x% j/+ Xs- ys~ Xs+ jys+ 

m-Es, A_B, T shapes 

Real D^' fractions 

Charge-flavor correlation 

Efficiency in the DP 

Background DP distributions 

_B~ -^ D'^^K^ cross-feed 

CP violation in Dn and BB 

Non-K* B- -> DK'^TT- decays 

Total experimental 0007 0004 0006 O004 O0T9 O009 OOTf O008 O029 O027 O029 O036" 



0.001 


0.001 


0.001 


0.001 


0.004 


0.006 


0.008 


0.004 


0.006 


0.003 


0.004 


0.002 


0.002 


0.001 


0.001 


0.001 


0.003 


0.003 


0.002 


0.002 


0.004 


0.001 


0.001 


0.001 


0.003 


0.003 


0.002 


0.001 


0.005 


0.005 


0.008 


0.002 


0.001 


0.001 


0.003 


0.001 


0.003 


0.001 


0.003 


0.001 


0.001 


0.001 


0.001 


0.001 


0.003 


0.001 


0.002 


0.001 


0.005 


0.002 


0.005 


0.003 


0.003 


0.002 


0.004 


0.004 


0.010 


0.004 


0.007 


0.002 


- 


- 


- 


- 


0.002 


0.003 


0.009 


0.002 


- 


- 


- 


- 


0.002 


0.001 


0.001 


0.001 


0.017 


0.001 


0.008 


0.004 


0.017 


0.002 


OOll 


0.001 


- 


- 


- 


- 


- 


- 


- 


- 


0.020 


0.026 


0.025 


0.036 



TABLE IL Summary of the main contributions to the Z)° decay amplitude model systematic uncertainty on the CP param- 
eters. We evaluate the different contributions using a similar, but not identical, procedure to that adopted in our previous 
analysis [9]. The reference D" decay amplitude models and parameters are used to generate 10 data-sized signal samples of 
pseudo-experiments of D*"*" — >• D" tt"*" and D*^ — > D ' ir^ events, and 10 B^ — >■ D^*' K^ and B^ — >■ D K*^ signal samples 100 
times larger than each measured signal yield in data, with D" -^ K^K^h^ . The CP parameters are generated with values in the 
range found in data. We then compare experiment-by-experiment the values of zC and Zszp obtained from the CP fits using the 
reference amplitude models and a set of alternative models obtained by repeating the D*' — >■ K^h^h^ amplitude analyses on the 
pseudo-experiments with alternative assumptions [12]. This technique, although it requires large computing resources, helps 
to reduce statistical contributions to the amplitude model uncertainties arising from changes in sensitivity between alternative 
models (e.g. alternative K-matrix solutions and P- vector mass dependence in the tttt S-wave parameterization). A variety of 
studies using data have been performed to test the consistency of the results using this procedure with those obtained in our 
previous analysis, where the alternative models were obtained by repeating the D^ — )• K%K^h~ amplitude analyses on data. 
Nevertheless, the largest decrease in the amplitude model uncertainty compared to our previous result is a consequence of 
the improvements in the experimental analysis of tagged D mesons [12] , which is reflected in smaller experimental systematic 
uncertainties on the D" decay amplitudes (variations of the reconstruction efficiency across the DP, modeling of the DP distri- 
butions for background events containing misreconstructed D mesons, mistag rates, etc.), and thus smaller amplitude model 
uncertainties on the CP parameters. 



Source 


X- 


y- 


x+ 


y+ 


x*_ 


y- 


xX 


y*+ 


Xs- 


Vs- 


Xs + 


ys+ 


Mass and width of Breit-Wigner's 


0.001 


0.001 


0.001 


0.002 


0.001 


0.002 


0.001 


0.002 


0.001 


0.002 


0.001 


0.002 


TTTT S-wave parameterization 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.002 


0.001 


0.001 


0.001 


0.002 


Kit S-wave parameterization 


0.001 


0.004 


0.003 


0.008 


0.001 


0.006 


0.002 


0.004 


0.003 


0.002 


0.003 


0.007 


Angular dependence 


0.001 


0.001 


0.002 


0.001 


0.001 


0.001 


0.001 


0.002 


0.002 


0.001 


0.002 


0.001 


Blatt-Weisskopf radius 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.002 


0.001 


Add/remove resonances 


0.001 


0.001 


0.001 


0.001 


0.001 


0.002 


0.001 


0.001 


0.001 


0.001 


0.001 


0.002 


DP efficiency 


0.003 


0.002 


0.003 


0.001 


0.001 


0.001 


0.001 


0.001 


0.004 


0.002 


0.003 


0.001 


Background DP shape 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


0.001 


Mistag rate 


0.003 


0.003 


0.002 


0.001 


0.001 


0.001 


0.001 


0.001 


0.003 


0.003 


0.001 


0.001 


Effect of mixing 


0.003 


0.001 


0.003 


0.001 


0.001 


0.001 


0.001 


0.001 


0.003 


0.001 


0.003 


0.001 


DP complex amplitudes 


0.001 


0.001 


0.001 


0.002 


0.001 


0.001 


0.001 


0.002 


0.002 


0.001 


0.001 


0.002 


Total D" decay amplitude model 


0.006 


0.006 


0.007 


0.009 


0.002 


0.007 


0.003 


0.006 


0.007 


0.006 


0.006 


0.008 



,(*) 



.(*). 



„(•) 



TABLE III: CP- violating complex parameters zi- = xl- +«J/^ and Zs^ — Xs^ +iya^ 



as obtained from the CP fit to KsTt^tv 
and KsK'^K~ final states separately. The first error is statistical, the second is the experimental systematic uncertainty and 
the third is the systematic uncertainty associated with the D" decay amplitude models. These results yield for the weak phase 



7 = 


(6I+i?)° {3,3}° and 7 =(87 


^37)° {8,3}°, respectively. 








^i!' 


T + TT- 


K'^K+K 


- 




Real part (%) 


Imaginary part (%) 


Real part (%) 


Imaginary part (%) 


z_ 


3.6 ± 4.6 ±0.9 ±0.6 


6.7 ± 4.9 ±0.4 ±0.6 


12.6 ± 7.6 ±1.5 ±0.5 


4.4 ± 11.7 ± 2.3 ±1.2 


z+ 


-8.3 ± 4.1 ±0.7 ±0.8 


-0.8± 4.9 ±0.4 ±1.0 


-19.0 ± 8.7 ±2.2 ±0.5 


- 2.0 ± 18.8 ± 6.0 ±1.5 


zl 


-8.9 ± 5.8 ±1.7 ±0.2 


-7.1± 6.8 ±0.9 ±0.7 


-17.0 ±11.0 ±2.0 ±0.4 


8.8 ± 17.1 ± 3.6 ±1.2 


4 


15.4 ± 5.9 ±1.4 ±0.4 


-3.6± 8.7 ±0.9 ±0.7 


11.7 ±11.7 ±4.2 ±0.4 


- 2.5 ± 16.4 ± 1.9 ±0.5 


Zs- 


12.8 ±10.5 ±3.4 ±0.8 


12.0 ±10.5 ±2.5 ±0.6 


-11.7 ±20.8 ±8.2 ±0.8 


14.3 ± 22.4 ± 9.5 ±1.4 


Zs+ 


-9.6 ± 9.2 ±3.2 ±0.8 


3.8 ±10.9 ±3.7 ±0.9 


-36.6 ±20.1 ±5.8 ±0.6 


-17.1 ± 39.9 ± 13.5 ± 1.7 



TABLE IV: Statistical correlation coefficients for the vector z of measurements, (in order) x_, j/-, x+, y+, x'L, yl, x\, y*^, 
Xs-, ys-, Xs+, ys+, as obtained from the CP fit to KsTv'^it^ and KsK'^ K~ final states (upper panel), and to Ksn^ir" (bottom 
left panel) and K'^K'^K^ (bottom right panel) separately. Only lower off-diagonal terms are written, in %. 
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-2 100 
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5 100 

100 

13 100 

100 

0-10 100 
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100 

15 100 

100 

0-19 100 

100 
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TABLE V: Experimental systematic correlation coefficients for the vector z of measurements, (in order) X-, j/_, x+, y+, x'L, 
j/1, x1, y+, Xs-, ys~, Xs+, ys+, defined as p^ — dj/^CaCjj, where dj — {7, — 7,hcst)i{'z — ^hast)j , with Zbest the vector of 
best measurements, for K'^-k'^-k~ and K^K^ K^ final states (upper panel), and K'^-k^-k^ (bottom left panel) and K%K'^ K~ 
(bottom right panel) separately. Only lower off-diagonal terms are written, in %. 



/lOO 
-23 100 
-2 18 100 



\ 



-5 
-32 
-15 
-12 

12 

-19 

3 

-20 

-1 



-5 
-10 
13 
-9 
17 
-8 
-6 

-6 



100 
■13 5 100 

-5 -2 100 
5 44 -60 100 
-9 -50 19 -44 100 
-7 49 25 -25 100 

-5 -2 5 63 100 
16 -14 73 65 100 



-2 
-9 
16 
■15 
-3 
-7 
-3 



/lOO 

13 100 

-3 -8 100 

-9 -6 2 100 

-28 -27 16 10 100 

3 2-5 -10 3 100 

-4 -1 10 5 16 -72 100 

16 7 -12 -16 -48 39 -57 100 

-11 -9 5 -6 40 -14 22 -29 100 



-3 -8 
33 
-3 



-6 -16 



2 64 100 



-11 -12 
-2 -6 



28 -13 10 -21 74 70 100 



3 2 66 97 74 100/ 
/lOO 

-72 100 

-78 59 100 

-91 78 88 100 

-41 30 17 34 100 

10 -11 -12 -5 -52 100 

-2 -11 -7 -1 60 -52 100 

-20 47 1 19 44 -42 46 100 

-13 5-8 3-1-2 3 9 100 



12 -6 6 



7 



-6 -98 100 



-25 32 40 36 7 15 -12 3 -80 81 100 



3 71 96 82 100/ \ 16 -21 -31 -26 -16 



-5 85 -85 -95 100 / 



TABLE VL Amplitude model systematic correlation coefficients for the vector z of measurements, (in order) X-, y_, x+, y+, 
x'L, yl, x*f_, y'f_, Xs~, ys~, Xb+, ya+, defined as previously, for K'^-k'^-k" and K^K'^ K~ final states (upper panel), and K'^tv^-k^ 
(bottom left panel) and K^K^ K~ (bottom right panel) separately. Only lower off-diagonal terms are written, in %. 
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FIG. 4: 1 — CL as a function of (left panel) rs, rg, and ura, and (right panel) Sb, Sg, and Ss, for B^ — s- DK^ , B^ -^ D* K^ , 
and B^ -^ DK*^ decays, including statistical and systematic uncertainties. The dashed (upper) and dotted (lower) horizontal 
lines correspond to the one- and two-standard deviation intervals, respectively. 




■0.5 0.5 



FIG. 5: (color online). Contours at 39.3% (dark) and 86.5% (light) 2-diniensional CL in the (a)(d) z^, (b)(e) z^, and (c)(f) 
Zsip planes, corresponding to one- and two-standard deviation regions (statistical only), for B~ (solid lines) and B'^ (dotted 
lines) decays, from the CP fit to the signal samples performed separately for (a)-(c) D — >■ A'sTt^tt" and (d)-(f) D — > K^K^ K~ 
decays. The results from the two subsets are statistically consistent. 
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FIG. 6: (color online). Contours at 39.3% (dark) and 86.5% (light) 2-diniensional CL in the (a)(c) z^,7r and (b)(d) 
planes, corresponding to one- and two-standard deviation regions (statistical only), for B~ (solid lines) and -B"*" (dotted lines) 
decays, from the CP fit to the B^ — > D^*'tv^ control samples performed separately for (a)(b) D — >■ KsTv'^tv~ and (c)(d) 



(*) 



0.01 



D — ^ KgK^ K~ decays. In this case we expect the z^,^ and z^^ contours close to the origin up to ~ 0.01, since r 
and the experimental resolutions are of the same order. Deviations from this pattern could be an indication that the DP 
distributions are not well described by the amplitude models [9]. The results from all the subsets are consistent with the 
expectations. Note the differences in scale when comparing to Fig. 5. 



